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未 来 30 年 川 东平 行 岭 谷 区 县 域 农田 SOC 动态 模拟 
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摘 要 : 以 位 于 川 东平 行 岭 谷 的 典型 县 垫 江上 县 为 研究 对 象 ， 探 讨 在 特定 气候 模式 下 未 来 30 a 研究 区 农田 土壤 有 机 碳 
(SOC) 及 其 动态 变化 ,为 研究 区 未 来 耕地 可 持续 利用 与 管理 提供 数据 支持 和 科学 依据 。 利 用 生物 地 球 化 学 模型 DNDC, 选取 
IPCC AR4 报告 中 的 BCCR_BCM 2.0 的 Bi 模式 ， 在 基于 研究 区 土壤 性 质 和 农业 管理 制度 等 建立 的 GIS 区 域 数据 库 的 支持 下 ， 
模拟 研究 区 2011 一 2041 年 SOC 动态 变化 。 结 果 表 明 : CDDNDC 模型 能 够 较 好 地 模拟 特定 气候 条 件 下 土壤 有 机 碳 及 其 动态 变 
化 ,模拟 值 和 观测 值 的 相关 系数 + 为 0.981, 达到 0.01 水 平 下 的 极 显著 相关 关系 ; 模拟 值 和 观测 值 的 RMSE 值 为 16%， 模拟 结 
果 较 好 。@) 未 来 30 a 研究 区 农田 0~20 cm 土 层 SOC 密度 和 储量 均 呈 显著 增加 态势 ， 单 位 面积 碳 增 量 2 637.07~8 091.55 kg(C) 
hm*， 增幅 为 10%~34%， 新 增 固 碳 量 2.7x105~8.3x105t 年 均 增 速 87.9~269.7 kg(C) hm*a!; @ 未 来 30 a 川 东平 行 区 县 域 农田 
土壤 总 体 呈 持续 碳 增 汇 状 态 , 研究 区 固 碳 、 丢 碳 以 及 相对 平衡 三 者 间 的 差异 将 逐渐 凸显 。 

关键 词 : 农田 土壤 ; 土壤 有 机 碳 ; DNDC 模型 ; 动态 模拟 ; 固 碳 量 ; 川 东平 行 岭 谷 区 
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Dynamic simulation of farmland SOC in parallel ridge valley area in Eastern 


Sichuan in the next 30 yearS 


ZHANG Lin, SHAO Jingan” 


(Department of Geography and Tourism, Chongqing Normal University / Key Laboratory of Surface Processes and Environmental Remote 


Sensing in the Three Gorges Reservoir Area, Chongqing 400047, China) 


Abstract: The relationship between greenhouse gas emission reduction and soil carbon sequestration has become the focus of global 
climate change research amidst the worsening global climate change. As a core interface of terrestrial ecosystem, farmland soils have a 
huge potential to sequester carbon. Exploring the response and feedback of the dynamic balance of farmland soil organic carbon (SOC) to 
global climate change is significant in understanding the relationship between global carbon cycle and agricultural ecosystems. Dianjiang 
County is a typical representation of parallel ridge valley area in Eastern Sichuan. It is an important grain and oil producing areas in 
Chongqing. In the process of transformation from traditional to modern agriculture, the area has been faced with double pressures of food 
security and increasing farmland SOC. Therefore in this region, it is very important to study farmland SOC and its dynamic change using 
specific climate model in the next 30 years. It provides data support and scientific basis for sustainable utilization and management of 
future cultivated lands in the study area. Supported by GIS regional database established for soil properties and agricultural management 
system, this study Simulated the dynamic change in SOC in the study area for the period 2011-2041 using the biogeochemical model 
DNDC. The study also used the Bi climate model within BCCR BCM 2.0 in the IPCC AR4 report. The results showed that: (1) The 
DNDC model fairly simulated farmland SOC and dynamic change under Specific climate conditions. The correlation coefficient between 
the simulated and measured values was 0.981, which was significant at the 0.01 confidence level. The RMSE value of the simulated and 


measured values was 16%, indicating better simulated results. (2) The SOC density and storage in the 0-20 cm soil layer in the study area 
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significantly increased for the next 30 years. The increase in carbon per unit area was 2 637.07-8 091.55 kg(C) hmz2 or 10%-34%, that in 
carbon sequestration Was 2.7x 10°5—8.3x107 t, and that in average increase rate was 87.9-269.7 kg(C) hm2.a-. (3) For the next 30 years, the 
farmland soil in parallel ridge valley area in Eastern Sichuan was generally under continuous carbon sink. In the study area, the differences 
among carbon sink, carbon loss and carbon relative balance gradually increased for the next 30 years. 

Keywords: Farmland soil; Soil organic carbon; DNDC model; Dynamic simulation; Carbon sequestration; Parallel ridge valley; Eastern 


Sichuan 


随 着 全 球 气候 变化 的 日 趋 严峻 , 温室 气体 减 排 和 土壤 碳 的 固 存 关 系 成 为 全 球 气候 变化 研究 的 热点 ! 
土壤 固 碳 是 一 项 重要 的 固 碳 减 排 措施 , 研究 土壤 有 机 碳 (SOC) 的 动态 平衡 对 全 球 0 
对 深入 了 解 全 球 碳 循环 和 农业 生态 系统 具有 重大 影响 。 累 计 文献 发 现 , 早期 的 研究 主要 是 根据 剖面 数据 进 
行 土壤 碳 储 量 的 估算 如 Post 由 根据 全 球 2 696 个 土壤 剖面 估计 全 球 土壤 有 机 碳 约 为 1.4.102t 土壤 碳 储 
量 约 占 陆地 生态 系统 碳 储量 的 2/3; 王 绍 强 等 站 依据 第 2 次 全 国土 壤 普 查 统计 资料 ,估算 中 国 陆地 总 SOC 储 
量 约 为 9.24:10%* g。 且 伴随 研究 的 进一步 深入 ， 越 来 越 多 的 学 者 利用 土壤 数据 库 、3S 技术 优化 碳 库 估算 和 
时 空 甄别 四 ,侧重 于 SOC 时 空 分 布 特征 与 影响 因素 方面 探索 。 揣 小 伟 等 四 借助 GIS 对 江苏 省 1985 一 2005 
年 表层 土壤 碳 密 度 进 行 核 算 , 探讨 了 土地 利用 变化 对 土壤 碳 储量 的 影响 ; 赵 永 存 等 四 运用 多 元 线性 回归 、 泛 
克 里 格 和 回归 克 里 格 3 种 方法 预测 了 河北 省 SOC 密度 的 空间 分 布 格局 。SOC 对 环境 条 件 的 动态 响应 和 自 
os 其 周转 过 程 非常 缓慢 ， 由 于 时 间 和 空间 的 局 限 , 现 有 的 田间 试验 无 法 揭示 不 同 

境 条 件 下 农田 SOC 的 动态 变化 。 近 年 来 模型 的 应 用 在 一 定 程度 上 解决 了 SOC 演变 的 模拟 与 预测 难题 。 
0 亮 等 名 利用 DNDC 模型 预测 了 新 疆 灰 漠 土 农田 SOC 的 动态 变化 ， 李 悦 等 0 借助 DAYCENT 模型 对 东 
北 地 区 不 同 农 作 管 理 措施 下 农田 SOC 未 来 变化 进行 了 模拟 研究 。 气 候 在 时 空 上 的 复杂 性 与 异 质 性 以 及 人 
类 活动 的 频繁 干扰 ,多 重 因素 使 得 SOC 与 全 球 气候 变化 的 相互 作用 在 现 有 研究 中 仍 存 在 较 大 的 争议 性 , 预 
六 未 来 气候 变化 下 SOC 动态 变化 对 深入 理解 陆地 生态 系统 与 全 球 气 候 变化 之 间 的 相互 作用 积极 重要 的 

意义 。 就 目前 而 言 模 型 模拟 的 途径 是 研究 未 来 土壤 碳 储 藏 变化 的 重要 手段 ，DNDC 
(DeNitrification- DeComposition) 模 型 是 利用 数学 模型 研究 生物 地 球 化 学 循环 过 程 ,模拟 陆地 生态 系统 中 碳 、 
所 和 水 循环 的 动态 模型 品 ， 现 已 被 广泛 应 用 于 农业 土壤 肥力 与 温室 气体 排放 的 预测 研究 中 ， 且 逐渐 被 认为 
是 目前 较为 典型 的 g 对 土壤 固 矶 强度 、 多 种 气体 排放 等 过 程 进 行 详尽 模拟 的 模型 。 

垫 江 县 是 川 东平 行 岭 谷 区 的 典型 代表 ,是 重庆 市 重要 粮油 产 区 之 一 , 在 传统 农业 向 现代 农业 转型 过 程 
1， 该 区 面临 担负 粮食 安全 和 增加 农田 SOC 的 双重 压力 。 鉴 于 此 , 本 文 以 垫 江 县 为 研究 区 , 在 课题 组 对 研 
究 区 1980 一 2011 年 SOC 动态 变化 号 研究 的 基础 上 , 根据 第 2 次 土壤 普查 和 2011 年 SOC 实测 值 对 DNDC 
模型 进行 参数 校正 , 选取 IPCC AR4 报告 中 的 BCCR_BCM 2.0 的 Bi 模式 ,开展 2011 一 2041 年 研究 区 SOC 
的 动态 变化 模拟 ， 结 果 有 利于 将 类 似 区 域 的 耕地 利用 提升 到 区 域 应 对 全 球 气候 变化 和 保障 粮食 安全 战略 的 
实施 层面 ,为 未 来 耕地 可 持续 利用 的 实践 提供 数据 支持 和 科学 依据 。 


1 材料 与 方法 
1.1 区 域 概况 


研究 区 域 垫 江 县 位 于 重庆 东北 部 川 东 平行 岭 谷 区 的 核心 区 (29°38'~30°31'N，107°13'~107°38 正 )， 面 积 
1518 km*( 图 1)。 地 处 华 莹 山脉 东部 ， 以 构造 地 貌 为 主 ,东西 部 山岭 务 峙 , 山 间 槽 地 交错 ， 中 部 高 滩 河 纵 贯 ， 
系列 溪 河 、 小 冲 沟 、 平 坝 镶 骨 , 形成 绵 豆 起 伏 的 丘陵 , 海拔 320~1 183 mm。 地 势 北 高 南 低 ， 以 平 坝 和 浅 丘 为 
] 气候 区 , 年 均 温 17 'C, 年 均 降 水 量 1 183 mm, 无 霜 期 289 d, 适宜 于 多 种 作物 种 植 。 
土壤 类 型 以 紫色 土 、 水 稻 土 、 黄 壤 和 冲积 土 为 主 , 水 稻 士 面积 4.36 万 hm”， 占 耕地 资源 的 42.7%。 研 究 区 
的 平坦 地 形 和 便捷 的 交通 , 使 得 该 区 成 为 重庆 市 传统 农业 同 现代 农业 转型 的 重点 县 , 为 配合 这 一 转型 ,该 
区 被 纳入 重庆 市 高 标准 基本 农田 建设 的 重点 县 , 因此 该 区 在 川 东 平行 岭 谷 区 具有 显著 的 代表 性 。 为 此 ， 该 
区 一 方面 要 担负 起 重庆 粮仓 的 重要 功能 , 通过 加 大 对 农业 的 基础 设施 建设 (如 土地 平整 、 路 沟 ( 渠 ) 池 由 等 ) 
和 农用 物资 (化 肥 、 有 机 肥 、 机 械 等 ) 的 投入 , 在 提高 劳动 生产 率 、 增 加 耕地 产 出 等 方面 获得 了 很 好 的 成 效 ; 另 
一 方面 这 一 投入 也 会 对 农田 土壤 产生 较 大 程度 的 扰动 , 如 土地 平整 、 机 械 翻 耕 等 .从 农田 SOC 的 保持 看 , 在 
这 一 区 域 具有 改善 与 扰动 的 双重 作用 ,是 研究 平行 岭 谷 区 传统 农业 向 现代 农业 转换 中 SOC 变化 的 代表 性 
区 域 。 


下 


201711.02086v1 


chinaXiv 


12 


Dem/m 


项 2750 
Low :16 


加 蕉 江 Dianjian 
0 4 8 


‘mer— KM 


ng County 


河流 River 
。 垫 江 采 样 点 
Soil sampling point 


ChinaXiv 合 作 期 刊 


了 花色 
有 黄 壤 土 Yellow soil 


图 1 研究 区 垫 江 县 位 置 、 数 字 高 程 模型 土壤 采样 点 和 耕地 土壤 类 型 分 布 
Fig.1 Location of Dianjiang County and its DEM, soil sampling site and farmland soil distribution 


数据 获取 与 设置 
1) 土 壤 数 据 : 1980 年 全 国 第 


2 次 土壤 普查 数据 来 源 于 垫 江 县 农业 局 ,提取 并 建立 土壤 空间 数据 库 所 需 


要 的 属性 数据 , 缺少 的 部 分 用 《 垫 江 县 土 种 志 》 补 缺 。2011 年 实测 样 点 分 析 数 据 ， 以 第 2 次 土壤 普查 原 采 


样 点 为 参考 ， 用 原始 采 柱 
样 点 展开 再 定位 , 在 原 采 样 点 及 其 附近 寻找 同 


提供 


根据 模型 需要 将 气象 数据 设置 成 晶 


人 NU 


采样 点 理化 性 状 的 最 终 值 。 


2) 气 象 数据 : 1980 一 2011 年 逐 


收 害 


1 和 施肥 等 种 植 制度 和 管理 措施 ,根据 执 洒 
4) 图 形 数据 : 1980 年 和 2011 年 的 耕地 资源 分 布 图 、 


点 记录 中 的 点 位 


述 (如 土壤 类 型 、 


地 名 、 河 流 、 
一 土 种 类 型 布设 土壤 采 检 


海拔 等 ), 在 ArcGIS 10.2 中 对 采 
点 。 依据 热 江 县 第 2 次 土壤 普查 数 


、 番 粒 、pH、 容 重 等 ), 将 计算 出 的 对 应 值 作 算术 平均 记 作为 


日 气象 数据 (降水 、 气 温 等 ) 由 重庆 市 气象 局 、 垫 江 县 及 周围 6 个 气象 站 点 
,2011 一 2041 年 的 逐日 气象 数据 采用 IPCC AR4 报告 中 的 BCCR_BCM 2.0 的 Bi 气候 模式 进行 设置 。 


最 高 气温 、 日 最 低 气温 及 日 平均 降水 量 。 
3) 作 物 数据 : 来 源 于 垫 江 县 统计 年 鉴 ， 选择 播种 面积 比较 大 的 6 种 种 植 制度 (冬小麦 - 夏 玉米 、 水 稳 -油菜 、 


[县 物候 记载 和 野外 调 


行政 区 划 图 


10.2 与 SigmaPlot 处 理 图 形 , Microsoft Excel 和 SPSS 17.0 统计 计算 。 


1.3 


建立 在 理论 ee 个 理 


DNDC 模型 参数 设置 


DNDC 模型 版 本 为 9.5。DNDC 模型 是 将 碳 、 


统 的 特征 。 然 而 ,尽管 模型 的 每 
但 是 土壤 碳 、4 ns zs 差异 怕 
施肥 等 ) 加 之 模型 原本 采用 的 参数 人 


验证 ， 


均 昌 


玉米 -红薯 、 玉 米 -黄豆 、 状 椒 -土豆 、 水 稻 - 豆 角 )， 进 行 模型 作物 子 模块 的 参数 设置 。 作 物 的 播种 、 灌 溉 、 


氮 循 环 的 物理 、 


需 通过 实测 对 部 分 影响 碳 、 


一 出 


行 修改 ， 如 水 稳 、 


些 重要 参数 特别 是 土壤 输 


De 


拟 


里 参数 ， 然 后 将 所 需 气 象 、 土 壤 及 田间 管理 
和 SOC 动态 变化 结果 与 田间 实测 
结果 与 田间 实测 值 拟 合 程度 达到 最 佳 。 


过 程 都 需 


以 到 


表 1 DNDC 模型 部 分 本 地 化 参数 值 


结果 的 拟 合 度 ,并且 还 需 多 次 1 


i 入 库 中 的 土壤 参数 进行 修改 和 本 地 化 ( 表 
小 麦 、 玉 米 、 大 豆 作 物 各 部 分 CON、 作 物 各 部 分 比例 等 。 首 先 设 置 模型 中 土壤 和 作物 生 
E 等 参数 输入 模型 ， 并 运行 ， 比较 模型 模拟 的 各 种 作物 类 型 生物 量 
则 整 DNDC 模型 中 的 部 分 参数 ， 以 使 模型 模 


查 数据 进行 综合 设置 。 
日 热 江 县 国土 局 提供 。 采 用 ArcGIS 


化 学 及 生物 过 程 用 计算 机 技术 组 织 起 来 ， 
想 化 的 生态 系统 ， 需 要 用 实测 数据 验证 其 是 否 具备 真实 生态 系 
论 或 观测 数据 为 依据 ,并 在 众多 农田 生态 系统 中 得 到 了 
FE， 而 且 模型 涉及 参数 繁多 (如 人 气象、 土壤 、 


决 省 值 主要 来 源 于 欧美 地 区 的 观测 结果 , 在 任 
氮 循 环 过 程 的 重要 参数 进行 本 地 化 及 模拟 结果 验证 。 为 此 , 研究 需要 对 模型 的 


1)， 对 模型 


何 特定 的 时 空 尺度 应 用 时 ， 


中 己 有 作物 的 默认 参数 进 


Table 1 Localized values of some parameters of DNDC model 


Parameter 
Clay fraction 
Wilting 
Porosity 


SOC decrease rate 


中 文 对 照 


Chinese reference 


黏 粒 含量 


凋 莹 点 
孔隙 度 
表层 以 下 土 


二 有 


Default value 


缺 省 值 本 地 化 值 依据 
Localized values Basis 
0.09 0.22 实测 Observed 
0.20 0.18 实测 Observed 
0.435 0.476 实测 Observed 
0.50 0.27 实测 Observed 
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below topsoil 机 质 减少 速度 


以 乡 (镇 ) 为 基本 模拟 单元 , 将 土壤 图 斑 数 据 库 与 热 江 县 行政 区 划 图 琶 加 , 产生 新 的 村 界 图 斑 数 据 库 。 按 
DNDC 区 域 模式 缺 省 数据 库 的 构建 方法 ， 以 图 斑 为 基本 模拟 单元 和 控制 区 域 , 用 地 统计 分 析 (Geo-statistical 
Analysb) 对 乡 (镇 ) 村 界 图 斑 的 土壤 属性 (初始 SOC、 秋 粒 含量 、pH、 容 重 等 ) 进 行 克 里 金 插 值 ， 提 取 土 壤 属性 
! 各 属性 值 的 均值 、 最 大 (小 ) 值 。 用 ArcGIS 10.2 空间 插值 功能 得 到 模拟 单元 内 各 作物 氮肥 和 农家 肥 投 入 及 
秸秆 还 田 的 空间 分 布 ， 形 成 县 域 土壤 数据 库 。 

DNDC 模型 与 GIS 数据 库 通 过 唯一 ID 识别 码 连接 ,用 垫 江 县 多 年 观测 数据 率 定 DNDC 模型 参数 , 开 

展 农田 SOC 由 点 位 到 区 域 的 模拟 , 并 在 点 位 和 区 域 尺 度 对 DNDC 模型 进行 验证 , 用 区 域 模拟 结果 与 实测 数 
据 建 立 环境 关系 模型 ， 对 残 差 进行 拟 合 ， 以 实现 SOC 模拟 尺度 的 扩展 。 
14 数据 分 析 
影响 SOC 密度 的 因素 涉及 地 形 、 上 土壤 理化 性 状 、 施 肥 、 管 理 、 生物 及 气候 等 0 
1)SOC 密度 初始 值 的 高 低 直接 影响 SOC 的 变化 , 其 值 越 大 ,变化 速率 越 慢 中 ,为 此 , 2011 年 的 SOC 密 
度 是 分 析 未 来 30 a 热 江 县 SOC 密度 动态 变化 的 首要 考虑 因素 。 
PN 犬 的 空间 分 异 有 很 大 影响 ,尤其 是 海拔 和 坡度 。 坡 度 大 小 和 海拔 高 低 决 
定 人 为 扰动 胁迫 的 程度 ,影响 土壤 成 士 过 程 ,进而 影响 SOC 的 空间 分 异 。 
3) 土 壤 理化 性 状 是 土壤 环境 中 较为 稳定 且 变化 较 慢 的 因素 ,在 影响 SOC 动态 变化 的 影响 因素 识别 过 程 
一 1， 是 必需 考虑 的 重要 方面 物理 性 状 中 如 容重 和 砾石 含量 对 SOC 的 影响 较 显 著 ; 化 学 性 状 中 除 pH A 有 
PP 效 磷 、 全 ( 氮 、 磷 和 钾 )、C/N 比 等 , 均 对 SOC 有 正 向 影响 上 ]。 这 类 影响 因子 中 , 主要 选择 容重 、 砾 石 含 
质地 、 有 效 磷 、 全 ( 氮 、 磷 和 钾 )、C/N 比 等 。 

4) 化 肥 ( 磷 、 氮 和 钾 ) 能 增加 土壤 养分 元 素 ， 有 机 肥 施 用 会 增加 有 机 质 含量 。 

5) 作 物产 量 高 意味 着 作物 生长 好 ,其 产生 的 枯 枝 落叶 较 多 ,会 增加 土壤 的 有 机 质 含量 ， 留 着 和 秸秆 还 
田 亦 是 如 此 , 直接 增加 土壤 有 机 质 ,进而 增加 SOC 含量 5。 

6) 气 象 数据 对 模型 输出 结果 有 重要 影响 , 如 降雨 量 、 大 气温 度 等 对 土壤 剖面 的 温度 、 湿 度 、 氧 化 还 原 
电位 、 碳 、 氧 等 有 直接 影响 。 研 究 将 未 来 30 a 的 土壤 管理 设置 为 常规 管理 模式 , 气象 因素 及 受气 象 因素 影 
自 的 作物 生长 状况 是 未 来 30 a 影响 SOC 储量 变化 的 主要 影响 因子 。 
以 点 位 模拟 为 基础 、 实 测 值 为 依据 率 定 模型 参数 ,开展 区 域 模拟 。 用 DNDC 模型 实施 点 位 和 区 域 模拟 ， 
须 事先 获取 区 域 所 有 格 点 所 需 输 入 数据 ,将 其 存储 于 特定 GIS 和 通用 数据 库 中 呈 。GIS 数据 涵盖 待 模拟 区 
所 有 多 边 形 和 格 点 信息 ， 如 地 理 位 置 、 气 象 数据 、 土 壤 性 质 、 种 植 制度 、 农 田 管理 等 。 通 用 数据 ， 如 植物 
生理 /物候 、 土壤 水 热 性 能 等 参数 。 为 支持 DNDC 模型 运转 , 首先 生成 以 村 为 最 小 单元 的 栅 格 , 建立 空间 数 
据 库 ， 并 假设 这 些 栅 格 性 质 是 均 质 的 ， 以 实现 对 各 个 栅 格 的 统一 模拟 ; 其 次 假定 每 个 栅 格 内 的 特征 是 一 致 
的 , 每 个 基本 单元 含有 且 仅 有 一 组 模型 所 需 输入 参数 , 再 在 每 一 基本 单元 模拟 的 基础 上 求 和 ”1。 在 县 域 
尺度 ， 首 先进 行 DNDC 点 位 模拟 ,以 此 率 定 模型 参数 ,再 以 垫 江 县 第 2 次 土壤 普查 、 自 然 环 境 等 作为 初始 
模拟 数据 ， 以 土壤 图 斑 法 , 结合 野外 实测 农田 管理 数据 建立 县 域 尺度 1 : 5 万 数据 库 ， 用 DNDC 模型 模拟 未 
来 30 a 县 域 尺度 SOC 情景 。 此 基础 上 , 用 热 江 县 40% 的 实测 样 点 作为 验证 数据 ,验证 DNDC 模型 模拟 结 
果 的 准确 性 。 
1.5 SOC 模拟 结果 与 精度 验证 

点 位 模拟 中 主要 选择 作物 生物 量 和 SOC 为 验证 指标 , 冬小麦 - 夏 玉 米 、 水 稳 -油菜 是 垫 江 县 典型 的 作物 
轮作 制度 对 其 进行 了 长 期 稳定 观测 和 土壤 理化 性 状 分 析 。 以 实测 土壤 理化 性 状 、 实 测 气象 资料 、 农 田 管 
理 信 息 为 DNDC 模型 的 输入 数据 , 将 模拟 得 到 的 产量 值 、SOC 与 实际 观测 值 对 比 验 证 模型 的 准确 性 。 据 实 
际 观 测 ， 以 冬小麦 常规 施肥 、 减 量 施肥 、 不 施肥 的 2003 年 、2008 年 、2011 年 的 实测 数据 及 水 稻 常规 施肥 
的 2005 年 、2009 年 、2011 年 的 实测 数据 对 模型 进行 验证 。 前 者 以 2003 年 为 初始 状态 , 用 冬小麦 - 夏 玉 米 旱 
地 SOC 实测 数据 ; 后 者 以 2005 年 为 初始 状态 ,用 水 稳 -油菜 水 田 SOC 实测 数据 ,耦合 气象 资料 和 农田 管理 
信息 运行 模型 。 

表 2 所 示 , 2008 年 小 麦 产量 模拟 值 与 实测 值 相对 误差 最 大 ， 为 -30.7%-~-12.2%, 均 相 对 误差 -16.2%， 模 
拟 结果 较 差 ， 其 中 不 施肥 处 理 相 对 误差 最 大 。2011 年 模拟 值 与 实测 值 较 接 近 , 介 于 -17.8%~-1.7%, 均 相 对 
误差 -2.7%， 模 拟 结 果 相 对 可 行 , 3 种 情景 下 模拟 值 丝 低 于 实测 值 。 表 3 可 知 , 2008 年 小 麦田 SOC 模拟 值 与 
实测 值 相对 误差 为 1.75%~16.75%, 均 相 对 误差 9.1%。2011 年 旱地 SOC 模拟 值 与 实测 值 相对 误差 为 
-2.18%-~9.66%， 均 相对 误差 3.1%。 模 型 模拟 初期 , 模拟 值 丝 低 于 实测 值 ， 原因 在 于 2003 年 前 垫 江 县 尚未 开 
展 试验 观测 ， 所 有 地 块 均 进 行 常规 施肥 处 理 , 之 后 才 实 施 差别 化 施肥 处 理 ， 从 而 使 得 2011 年 模拟 值 越 来 越 


Ei 


x 


接近 实测 值 


四 


且 伴 随 模 拟 时 间 的 延长 ,模型 越 趋 于 稳定 。 表 2 与 表 4 显示, 2009 和 2011 年 水 稳 模 拟 值 与 


党 


测 值 相对 误差 较 小 , 分 别 为 7.5%、-5.9%, 模拟 结果 较 理 想 。 2009 和 2011 年 水 田 SOC 模拟 结果 也 较 好 ,， 模 


拟 值 


与 实测 值 相对 误差 分 别 为 7.3% 和 -7.1%6。 


表 2 DNDC 模型 模拟 的 不 同 施肥 处 理 下 小 麦 产量 及 常规 施肥 下 的 水 称 产量 与 实测 产量 的 对 比 


Table 2 Comparison between Simulated and measured wheat yield under different fertilization treatments and rice yield unde 


conventional fertilization treatment 


be 年 份 ”施肥 处 理 i ee yield a 相对 误差 

Year Fertilization treatment (kg hm?) (kg: hm’) Relative error (%) 

小 考 。 2008 常规 施肥 Conventional fertilization 1 520 1334 _122 

Ws 减 量 施 肥 Balanced fertilization 1 564 1331 -14.9 

不 施肥 None fertilization 547 379 -30.7 

2011 常规 施肥 Conventional fertilization 1 616 1 557 _3.6 

减 量 施肥 Balanced fertilization 1 452 1 428 -1.7 

不 施肥 None fertilization 670 651 1 了 7 海 

水 稻 2009 常规 施肥 Conventional fertilization 3 032 2 804 _7.5 

we 2011 常规 施肥 Conventional fertilization 3 169 2 981 -5.9 


表 3 DNDC 模型 模拟 与 实测 的 2008 年 和 2011 年 小 麦田 不 同 施肥 处 理 土壤 有 机 碳 (SOC) 对 比 
Table 3 ”Comparison between Simulated and measured soil organic carbon (SOC) contents under different fertilization treatments in 
wheat field in 2008 and 2011 


2003 年 初始 


施肥 处 理 SOC 实测 值 
ertilizat 1 人 t Initial value Observed 
ertilization treatmen 0 i 


(gkg') (gkg’) 


常规 施肥 Conventional fertilization 4.78 5.14 
减 量 施肥 Balanced fertilization 4.52 5.73 
不 施肥 None fertilization 4.06 4.91 


2008 2011 
异 拟 值 、 实测 值 费 拟 值 、 
模拟 值 相对 误差 实测 值 模拟 值 相对 误差 
Simulated : Observed Simulated 
Relative Relative 
value ELLY value value LioOF 
1 
(gkg) %) (gkg) (gkg) (0%) 
5.23 1,73 5.51 5.39 —2.18 
6.69 16.75 5.80 6.36 9.66 
5.30 7.90 $5.02 5.07 0.10 


表 4 DNDC 模型 模拟 与 实测 的 2009 年 和 2011 年 水 稻田 常规 施肥 土壤 的 有 机 碳 (SOC) 含 量 对 比 
Table 4 Comparison between simulated and measured soil organic carbon (SOC) contents under conventional fertilization 
treatment in the paddy field in 2009 and 2011 


2005 初始 2009 2011 
SOC 实测 值 模拟 值 ee 实测 值 模拟 值 1 

Initial value Observed Simulated 相对 误 关 Observed Simulated 相对 误 关 
in 2005 Relative Relative 
1 value value error (%) value value error (%) 

本 二 0 0 

(gkg ) (Cgkg) (gkg) (gkg') (gkg ) 
4.9 6.23 6.7 7.5 8.64 8.02 -7.1 


利用 校正 后 的 土壤 性 质 和 作物 系数 ,模拟 执 江 县 1980 一 2011 年 农田 SOC, 将 模拟 均值 与 各 采样 点 观测 


性 甘 


一 2 xl 二、 


值 对 比 检验 ,本 研究 采 
1IRMSE 的 计算 公式 如 下 : 


相关 系数 (r) 和 相对 RMSE(root mean square error) 验 证 DNDC 模型 的 区 域 尺度 适宜 


100 


RMSE, “pps > ON SM,)/n (1D) 


式 中 : OBS; 为 观测 值 , SM; 为 模拟 值 ，n 为 样本 容量 


一 般 参 


考 为 : RMSE<10%, 表明 模拟 值 与 预测 值 


.RMSE 的 值 越 小 ,表示 模拟 值 与 实际 值 的 一 致 性 越 好 。 
到 性 


非常 好 ; 在 10%~20% 之 间 为 较 好 ; 在 20%~30% 之 间 


表明 模拟 效 


果 一 般 ; >30% 则 表明 模拟 值 与 实际 值 偏差 较 大 ,模拟 效果 较 差 ”。 


研究 区 农田 表层 SOC 含量 的 DNDC 模拟 值 与 田间 实测 值 表现 出 较为 一 致 的 变化 趋势 ， 二 者 相关 系数 7 
为 0.981, 达到 了 0.01 水 平 下 的 极 显 著 相 关 关 系 , 模拟 值 和 观测 值 的 RMSE 值 为 16%。 由 此 而 言 , DNDC 模 


型 模拟 的 效 


果 良 好 , 利用 DNDC 模型 模拟 热 江 县 2011 一 2041 年 农田 SOC 趋势 变化 是 可 行 的 。 
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2 结果 与 分 析 
2.1 SOC 密度 和 储量 动态 变化 

SOC 的 转化 是 一 个 长 期 而 又 复杂 的 过 程 ， 本文 利用 DNDC 模型 ,选择 IPCC AR4 报 告 中 的 BCCR_BCM 
2.0 的 Bl 模式 ,对 热 江 县 2011 一 2041 年 的 未 来 30 a 土壤 有 机 碳 进 行 模拟 预测 ， 模 拟 中 假设 模拟 时 间 内 耕作 
管理 方式 完全 相同 。 模 拟 结果 (图 2) 表 明 : 2011 一 2041 年 热 江 县 农田 0~20 cm 土 层 SOC 密度 和 储量 均 呈 显著 
增加 趋势 ,模拟 的 2041 年 的 SOC 密度 为 2011 年 原始 密度 的 110%~134%, SOC 密度 变化 为 2 637.07~8 091.55 
kg(C):hm”*。 不 同 土壤 类 型 SOC 增 减 幅度 有 所 差异 ,增幅 最 大 的 是 热 江 县 面积 第 4 的 冲积 土 , 增幅 达 
130%~151%。 增 幅 最 小 的 是 重庆 普遍 发 育 的 紫色 土 , 增幅 为 -9%~-7%， 面积 第 2 的 水 稻 士 和 面积 第 3 的 黄 
壤土 增幅 分 别 为 30%~42% 和 20%~ 30%。 

2011 一 2041 年 热 江 县 SOC 储量 变化 整体 上 与 SOC 密度 的 变化 趋势 一 致 ， 模 拟 2041 年 SOC 储量 约 为 
2011 年 原始 值 的 110%~168%， 新 增 固 碳 量 在 2.7x10”~8.3x10t 之 间 。 由 于 面积 的 影响 30 a 农田 不 同 土壤 类 
型 SOC 储量 变化 与 SOC 密度 变化 有 所 不 同 。 单 位 面积 碳 增 量 最 大 的 冲积 土 , 因 其 面积 比较 小 致使 新 增 固 

类 量 亦 较 小 ,而 占 耕 地 面积 最 大 的 紫色 土 和 第 2 的 水 稻 士 , 总体 决定 垫 江 县 农田 SOC 的 盔 亏 。 


SOC 密度 变化 a SOC 储量 变化 有 
可 23 Ne 6.Sx10: 1.9x10' 
SOC density change [kg(C)-hm’3] | >” 0S SOC storage change (0) 6 


O10 146 ~ 15 900 E23 108 — 31 io ON -L7xio ~ 1.4x10; mm -io ~— 4.1x10° 
5 5 0 4 
M1; "oo ~ 19 598 名 和 1.4x10; ~ 6.5x10 的 


图 2 2011 一 2041 年 农田 土壤 有 机 碳 (SOC) 密 度 (a) 和 储量 (b) 变 化 的 空间 分 布 格局 
Fig. 2 Spatial patterns of soil organic matter (SOC) density (a) and storage (b) during 2011-2041 


垫 江 县 未 来 30a 农田 SOC 密度 变化 总 体 呈 西北 、 东 南 高 ,中 间 低 的 空间 分 布 格局 ， 变 化 的 高 值 表现 在 
研究 区 西北 和 东南 方 , 变化 的 低 值 表现 在 研究 区 中 部 。 这 与 热 江 县 的 骨架 地 势 相 一 致 ,增幅 最 大 的 处 于 三 
大 山系 (明月 山 、 金 华山 和 黄 草山 ) 及 其 周围 区 ,其 次 是 中 部 高 滩 河 左上 游 流 域 、 大 沙河 右 下 游 流域 、 南 部 
及 稚 大 台地 区 ,最 低 的 是 明月 和 金华 两 山 所 夹 槽 谷 区 的 交错 地 带 。 垫 江 县 SOC 储量 变化 的 空间 格局 与 密度 
变化 的 分 布 格局 有 细小 差别 , 储量 增加 最 明显 区 位 于 北部 回 龙 河上 流 沿 岸 的 河谷 平 坝 区 和 金华 山西 南 段 ， 
其 次 是 三 大 山系 及 其 周围 区 , 增 量 最 小 的 是 明月 山 与 回 龙 河 所 夹 槽 谷 及 金华 山 与 高 滩 河 所 夹 槽 谷 区 。 

2.2 ”SOC 密度 和 储量 年 均 变 化 速率 及 丢 固 碳 分 布 特征 

未 来 30 a 垫 江 县 SOC 密度 和 储量 的 变化 速率 整体 均 呈 增加 趋势 (图 3)。2011 一 2041 年 农田 0~20 cm 土 
层 SOC 密度 模拟 年 均 增加 速率 为 87.90-269.72 kg(C)-hm*…a-!', 高 于 19892011 年 的 相应 年 均 增 长 率 72.11 
kg(C)-hm*a', 增幅 最 小 为 21.90%, 总 体 呈 较 明 显 的 碳 增 汇 效应 。 垫 江 县 农田 SOC 储量 变化 速率 与 相应 时 
间 段 SOC 密度 变化 速率 趋势 相 一 致 其 年 均 增长 速率 8.9x10”~2.8x104ta!, 远 高 于 1989_2011 年 的 相应 年 
均 增长 速率 7373.31 ta-"。 未 来 30 a 热 江 县 农田 SOC 密度 和 储量 的 变化 速率 在 不 同 土壤 类 型 间 差异 明显 。 
SOC 密度 年 均 增 长 速率 最 大 的 是 冲积 土 , 达 377.9~494.2 kg(C)-hm-*a-'; 其 次 是 水 稻 土 , 达到 288.0~401.2 
kg(C)-hm-*a-!; 再 次 是 与 水 稻 土 相差 不 大 的 黄 壤 土 , 达 260.0~491.2 kg(C)-hm-*…a-!; 最 小 的 是 紫色 土 , 增长 速 
率 为 -186.1~-51.9 kg(C)-hm*a-'。 与 1989 一 2011 年 农田 SOC 密度 在 不 同 土壤 类 型 间 的 年 均 增长 速率 [冲积 
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土 535.00 kg(C):hm-*a'、 黄 壤土 133.46 kg(C)-hm*a'、 


ChinaXiv 合 作 期 刊 


水 稳 土 218.48 kg(C)-hm*a 和 紫色 土 -50.64 


kg(C)-hm*a-] 比 较 , 冲积 土 增长 速率 有 所 降低 ， 黄 壤土 增长 速率 大 幅度 提高 ,水 稻 土 增长 速率 微弱 提高 ， 


紫色 土 增长 速率 出 现 降 低 。 


究 其 原因 ,土壤 所 能 容纳 SOC 量 是 有 限 的 ,过 去 30 a 冲积 土 SOC 密度 增长 速率 原本 较 大 ,达到 某 一 峰 


值 后 会 逐渐 减低 。 而 黄 壤 土 增长 速率 较 小 且 处 于 上 升 阶 段 ,换言之 , 3 


其 所 能 容纳 SOC 能 力 的 空间 较 大 。 水 


稻 土 因 人 为 长 期 耕作 扰动 , SOC 容纳 能 力 处 于 相对 稳定 状态 , 波动 较 小 。 紫 色 土 由 其 本 身 理化 性 状 决定 的 
SOC 容纳 能 力 难以 突破 。SOC 储量 变化 速率 也 呈 较 明显 地 增长 趋势 , 且 不 同 土壤 类 型 间 的 差异 既 受 SOC 
密度 变化 的 影响 , 也 受 对 应 土壤 类 型 面积 的 控制 。 与 水 稻 土 相 比 ,冲积 土 虽然 SOC 密度 年 均 变 化 速率 最 大 ， 
但 其 面积 较 小 ( 仅 为 水 称 土 的 0.23%), SOC 储量 变化 速率 仅 为 水 稻 土 的 SOC 储量 变化 速率 的 4%, 两 者 相差 
其 远 。 在 SOC 的 管理 中 , 不 但 要 考虑 SOC 密度 的 提升 空间 和 潜力 ,也 要 考虑 土壤 类 型 的 面积 占 比 ,两 者 结 


合 才 能 从 总 体 上 提升 SOC 含量 。 


以 2011 一 2041 年 SOC 储量 变化 量 占 2011 年 储量 比重 定义 : 固 碳 (>5%)、 技 碳 (<-5%)、 介 于 两 者 间 视 
为 相对 平衡 。 以 此 划分 ， 热 江 县 农田 土壤 固 碳 、 丢 碳 和 相对 平衡 面积 比分 别 为 58.14%、31.62% 和 10.24%， 从 
整体 来 看 (图 3), 未 来 30 a 垫 江 县 农田 土壤 固 碳 与 相对 平衡 的 面积 之 和 占 总 面积 的 203， 丢 碳 面积 仅 占 1/3， 
与 农田 SOC 储量 显著 增加 的 趋势 相 一 致 。 进 一 步 细 分 , 固 碳 幅度 为 30.19% 的 水 稻 士 ， 其 固 碳 面积 比较 大 ， 


为 84.60%, 面积 上 的 优势 使 得 其 成 为 垫 江 县 农田 碳 增 汇 关 键 。 黄 壤土 的 固 碳 面积 略 高 于 丢 砚 玫 


职 ， 占 比分 


别 为 57.00% 和 43.00%。 紫 色 土 固 碳 面 积 为 26.08%， 技 碳 面 积 为 58.72%， 相 对 平衡 面积 仅 15.20%， 因 其 丢 
碳 面积 显著 大 于 固 碳 面积 ， 这 在 很 大 程度 上 影响 县 域 尺 度 的 农田 碳 增 汇 效 应 。 垫 江 县 未 来 30 a 农田 土壤 的 


兢 源 / 汇 主要 由 水 稻 土 和 紫色 土 间 的 碳 效应 所 决定 ， 


且 对 所 有 土壤 类 型 来 讲 , 未 来 30 a 农田 SOC 储量 变化 


没有 表现 出 单一 的 丢 碳 、 回 碳 或 相对 平衡 格局 ， 而 是 丢 碳 与 固 碳 同时 并 存 或 丢 碳 、 固 碳 与 平衡 三 者 同时 存 


在 。 原 因 是 影响 农田 SOC 储量 变化 的 因素 繁多 , 同一 演化 时 期 内 ， 


种 因素 共同 作用 的 结果 。 


SOC 年 均 变 化 速率 
Annual average rate 
of SOC density |kg(C)-hm2a 1 1 123~2115 


标量 -7s ~ -209 加 于 215 ~ 448; 
0 4 8 
Km 


Em -209 ~ 1 123 


丢 ( 固 ) 碳 分 布 特征 
Distribution charater of 
carbon loss a 

a Carbon loss 
OM Fk carbon sink 


ater 0 
loss and carbon si 


nk 国有 千 衔 Carbon relative 
b: 


alance 
4 8 


Km 


图 3 2011 一 2041 年 土壤 有 机 碳 (SOC) 密 度 (a) 和 储量 (b) 的 年 均 变化 速率 空间 格局 
Fig.3 Spatial patterns of annual average rate of soil organic carbon (SOC) density (a) and carbon sink, loss and 
relative balance (b) during 2011-2041 


未 来 30 a 垫 江 县 农田 土壤 固 碳 主要 表现 在 研究 区 的 东南 、 西 北 前 
区 的 中 部 , 呈现 出 东南 、 西 北部 固 碳 和 中 部 丢 碳 的 分 布 格局 (图 3)。 固 碳 区 集中 于 和 锥 大 台地 中 南部 、 三 大 主 


主要 发 生 在 明月 山 与 下 缘 模 谷 的 交错 区 及 金华 山南 段 、 稚 大 台地 北 段 。 丢 碳 
坝 交 错 带 , 坡度 起 伏 较 大 , 不 利于 无 机 、 有 机 肥 的 施用 与 耕作 管理 ,作物 长 势 与 产量 相对 较 低 , 地 上 凋零 物 


固 碳 、 丢 碳 和 相对 平衡 状态 共存 是 由 多 


%， 丢 碳 和 相对 平衡 情况 主要 表现 研究 


体 山 系 ( 明 月 山 、 金 华山 和 黄 草山 ) 及 零星 散布 在 回 龙 河 所 在 槽 谷 平 坝 区 。 相 对 平衡 区 集中 在 西部 明月 山 右 
下 缘 槽 谷 区 、 高 滩 河 与 回 龙 河 所 夹 平 坝 区 、 高 滩 河 下 游 以 西 以 北 区 、 黄 草山 紧 靠 下 缘 的 以 西 以 北 区 。 


丢 碳 


区 普 融 属于 山系 下 缘 与 相 


三 | 


谷 平 


和 地 下 部 分 相对 较 少 ,致使 农田 SOC 储量 的 变化 整体 呈 丢 碳 状态 。 相 反 地 , 固 碳 区 为 平 坝 区 和 主体 山系 下 
缘 ， 前 者 利于 无 机 、 有 机 肥 投 入 及 耕作 的 管理 ,后 者 地 上 凋零 物 和 地 下 部 分 较 丰 富 , 则 利于 SOC 储量 的 提 
升 。 年 均 SOC 密度 变化 速率 在 空间 格局 上 呈 西 部 、 西 北部 高 于 东部 、 东 南部 , 更 高 于 东北 部 、 西 南部 (图 
3)。 大 体 分 为 三 级 : 西部 明月 山 和 东南 部 黄 草 山 是 变化 高 值 区 ; 高 滩 河 西北 、 大 沙河 东南 部 、 大 沙河 中 下 游 
河 段 区 是 变化 较 显 著 的 中 值 区 ; 县 城 以 东 ， 大 沙河 以 北 , 高 滩 河 与 箱 大 台地 所 夹 槽 谷 区 是 变化 低 值 区 。 


3 讨论 与 结论 
3.1 讨论 


影响 SOC 变化 的 因子 较为 复杂 ， 如 土壤 温度 、 湿 度 、 氧 化 还 原 、 水 土 流失 和 植物 生物 量 等 ,而 这 些 条 
件 又 受制 于 气息、 土质 、 植 被 、 人 类 活动 等 驱动 力 ， 是 影响 因子 在 组 合 上 的 任何 变化 都 会 改变 农田 SOC 的 
质 和 量 。 其 中 SOC 密度 初始 值 、 土 壤 全 所 密度 以 及 C/N 和 土壤 颗粒 组 成 对 土壤 有 机 碳 影响 较 大 。 程 先 福 
等 四 的 研究 发 现 SOC 密度 初始 值 对 农田 表土 SOC 密度 的 影响 最 大 且 呈 负 相 关 关 系 。 本 研究 也 表明 ，SOC 
密度 的 初始 值 [ 越 高 , 未 来 30 a 热 江 县 SOC 密度 年 均 变 化 速率 越 小 。 农 田 SOC 存在 容纳 上 限 ,，SOC 含量 越 
高 ,短期 内 越 难 显著 地 提升 ; 反之 ,通过 合理 耕作 、 施 肥 、 管 理 等 措施 , 在 短期 内 能 明显 提升 SOC 含量 请 和 。 
la SOC 密度 初始 值 高 的 北部 和 东北 部 , SOC 密度 的 提升 速率 低 于 相对 
较 低 的 中 部 、 南 部 和 东南 部 ,土壤 全 氮 密 度 及 C/N 对 农田 表层 SOC 密度 的 影响 仅 次 于 SOC 密度 初始 值 ， 与 
土壤 SOC 呈正 相关 关系 。 土 壤 全 氮 密 度 和 C/N 值 越 大 , 未 来 30 a 垫 江 县 SOC 密度 的 变化 速率 越 大 。 课 题 
组 通过 对 前 30 a 垫 江 县 SOC 动态 研究 也 发 现 全 氮 密 度 与 C/N 存在 一 定 依赖 性 ,土壤 中 N 的 积累 往往 伴随 
着 碳 的 增加 。 

国家 发 展 进 程 的 不 断 推 进 , 研究 区 耕地 面积 会 随 着 区 域 发 展 需求 发 生 相 应 的 改变 。 耕 地 面积 的 增加 或 
减少 ， 人 为 扰动 加 大 与 减少 都 将 导致 研究 区 土壤 SOC 密度 、 储 量 发 生变 动 放 5。 这 与 本 文 在 假设 研究 区 耕 
地 面积 不 变 情况 下 得 出 的 结论 有 出 入 。 运 用 模型 方法 对 农田 土壤 碳 库 及 固 碳 潜力 模拟 已 成 为 目前 研究 农田 
SOC 的 一 种 重要 方法 与 途径 ,本 文 利 用 生物 地 球 化 学 模型 DNDC, 仅 选 取 IPCC AR4 报告 中 的 BCCR_ 
BCM 2.0 的 Bl 模式 下 未 来 气温 和 降水 预测 结果 ,对 川 东平 行 岭 谷 区 垫 江 县 未 来 30 a 的 农田 SOC 的 变化 趋 
势 进行 模拟 , 略 显 单薄 。 要 准确 、 有 效 地 预测 本 研究 区 农田 SOC 的 变化 趋势 , 不仅 要 收集 更 多 的 样 点 以 及 
不 同 作物 的 土壤 数据 进行 模型 的 参数 率 定 与 验证 , 还 需 进行 更 多 预测 气候 模式 下 的 模拟 研究 。 此 外 ， 本 研 
究 是 在 以 统计 资料 与 实测 数据 为 基础 的 数据 库 的 支持 下 展开 的 ,其 结果 与 统计 资料 的 准确 性 及 实测 数据 的 
精度 有 密切 关系 , 这 也 会 导致 模拟 结果 与 实际 存在 一 定 偏 差 。 


3.2 结论 


pb 本 研究 采用 相关 系数 x 和 相对 RMSE(root mean square error)2 个 指标 验证 了 DNDC 模型 对 土壤 有 机 碳 

:= 的 模拟 效果 ,二 者 相关 系数 + 为 0.981, 达到 了 0.01 水 平 下 的 极 显著 相关 关系 ; 模拟 值 和 观测 值 的 RMSE 值 
为 16%, 说 明 模 拟 的 效果 较 好 。DNDC 模型 适合 于 大 尺度 、 长 时 间 范 围 的 农田 SOC 的 时 空格 局 模拟 ,伴随 
和 时 间 的 推移 ,模型 模拟 结果 趋 于 稳定 ， 且 与 实测 值 越 接 近 , 模拟 结果 越 理 想 。 
本 文 借助 DNDC 模型 对 川 东 平行 岭 谷 典型 县 未 来 30 a 农田 SOC 密度 进行 模拟 ,结果 显示 未 来 30 a 热 
江 县 农田 0~20 cm 土屋 SOC 密度 以 及 储量 都 将 呈现 出 明显 增加 的 态势 。 其 中 SOC 密度 由 2011 年 的 
23 798.68 kg(C)-hm 增 加 到 2041 年 的 26 435.70~31 890.23 kg(C)-hm2， 单 位 面积 碳 增 量 2 637.07~8 091.55 
kg(C)-hm”, 新 增 固 碳 量 2.7x105 t~8.3x10” t。 

未 来 30 a 研究 区 SOC( 密 度 和 储量 ) 变 化 速率 整体 上 均 呈 现 出 逐渐 增加 趋势。 DNDC 模型 模拟 出 
2011 一 2041 年 农田 0-20 cm 土 层 SOC 密度 年 均 增 速率 是 87.90-269.72 kg(C):hm*a， 远 远 高 于 1989 一 2011 
年 农田 0~20 cm 土 层 SOC 密度 年 均 增 长 值 72.11 kg(C)-hm*…a-',， 增长 幅度 高 达 21.90%, 总 体 上 将 呈现 出 比 
较 明 显 的 碳 增 汇 效应 

DNDC 模型 模拟 结果 显示 , 研究 区 未 来 30 a 丢 碳 面积 比例 为 31.62%、 固 碳 面积 比例 为 58.14%、 相 对 
平衡 的 面积 比例 为 10.24%。 固 碳 面积 大 于 丢 碳 面积 与 相对 平衡 面积 之 和 , 将 打破 课题 组 之 前 1980 一 2011 
年 热 江 县 固 碳 面积 与 竺 碳 、 平 衡 面 积 之 和 的 持平 状态 , 三 者 之 间 的 显著 性 将 逐渐 凸显 。 
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